A Derivation of the Classical Einstein-Dirac-Maxwell Equations 
From a Model of an Elastic Medium 



John Baker 
2221 Parnassus Ct; Hayward, CA 94542 

Abstract 

Starting from a model of an elastic medium, partial differential equations with the form of 
the coupled Einstein-Dirac-Maxwell equations are derived. The form of these equations describes 
particles with mass and spin coupled to electromagnetic and gravitational type of interactions. A 
two dimensional version of these equations is obtained by starting with a model in three dimensions 
and deriving equations for the dynamics of the lowest fourier modes assuming one dimension to be 
periodic. Generalizations to higher dimensions are discussed. 
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I. INTRODUCTION 



Dirac's equation describes the behavior of particles with mass and spin and how they 
couple to the electromagnetic field. The usual form of Dirac's equation (with Ti, and c set to 
unity) isi 

(7^'(9/, + m)^(x) = 

where \l/ is a spinor field, m is the particle mass and the gamma matrices satisfy 

The electromagnetic field is introduced by the minimal coupling prescription? — > D^, 
with 

where is the electromagnetic vector potential. Dirac's equation can be further coupled 
to gravity (at the classical level) and the equation then takes the formi^ii^i^'^i^ 

7'^[9^-r^ + ^eA^]^(x)+m$(x) =0 (1) 

where T ^ is known as the spin connection, is the electromagnetic vector potential, and m 
is the mass. The gravitational coupling enters through the modified dirac matrices t''^ which 
satisfy the anticommutation relation 



{7^r} = 2/^7'^^ (2) 
and the spin connection satisfies the additional constraint^i^ 



^ + ^f^V^ - V^i^ + i^V^ = (3) 

where Fj^j, are the usual Christoffel symbols. 

The above form of Dirac's equation describes the dynamics of the spinor field \1/ when 
coupled to the scalar fields and gravity. There are two additional equations which describe 
the dynamics of A^ and g^^, these are the Einstein field equation 



R^J.u — -R — T^u (4) 
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and Maxwell's equations 



y^^M^ = 4716^7^^^ (5) 

The Equations ([1]), (j4]) and ([5]) are collectively known as the Einstein-Dirac-Maxwell 
equations^i^. This paper will show that there exists a "Dirac-like" equation, with the 
form of Equation ([T]), that determines the dynamics of the lowest Fourier modes of an elastic 
medium and furthermore the dimensional reduction involved in the Fourier transform implies 
Equations and 

II. ELASTICITY THEORY 

The theory of elasticity is usually concerned with the infinitesimal deformations of an 
elastic body^'^»^'^»^. We assume that the material points of a body are continuous and can 
be assigned a unique label a. For definiteness the elastic body can be taken to be a three 
dimensional object so each point of the body may be labeled with three coordinate numbers 
a' with i = 1,2,3. 

If this three dimensional elastic body is placed in a large ambient three dimensional space 
then the material coordinates a* can be described by their positions in the 3-D fixed space 
coordinates with i = 1,2,3. In this description the material points a*(x^, x^, x^) are 
functions of x. A deformation of the elastic body results in infinitesimal displacements of 
these material points. If before deformation, a material point a° is located at fixed space 
coordinates x°^,x°^,x°^ then after deformation it will be located at some other coordinate 
. The deformation of the medium is characterized at each point by the displacement 

vector 

=x^ - x^' 

which measures the displacement of each point in the body after deformation. 

It is our aim to take this model of an elastic medium and derive from it equations of 
motion that have the same form as Dirac's equation. In doing so we have to distinguish 
between the intrinsic coordinates of the medium which we will call "internal" coordinates 
and the fixed space coordinates which facilitates our derivation of the equations of motion. 
In the undeformed state we may take the external coordinates to coincide with the material 
coordinates a* = x°*. The approach that we will use in this paper is to derive equations of 
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motion using the fixed space coordinates and then translate this to the internal coordinates 
of our space. 

We first consider the effect of a deformation on the measurement of distance. After 
the elastic body is deformed, the distances between its points changes as measured with 
the fixed space coordinates. If two points which are very close together are separated by 
a radius vector dx^^ before deformation, these same two points are separated by a vector 
dx'^ = + du^ afterwards. The square distance between the points before deformation is 
then ds'^ = (dx^^)"^ + {dx^'^Y + {dx^^y. Since these coincide with the material points in the 
undeformed state, this can be written ds'^ = {da^Y + {do?Y + {da^Y' . The squared distance 
after deformation can be written^^ ds'"^ = {dx^y+{dx'^y+{dx^y = J^iidx'^Y = J2i{da^+du^)'^. 
The differential element du^ can be written as du^ = Y,i ^da', which gives for the distance 
between the points 



ds- = E {da^ + E ^da'^) {da^ + E ^da'^ 

da" da" + -—da'' da" + V -—da' da' + V V -—r—rda''da' 
Y da'' Y ^« ki ^« ^« y 

= ^{5ik + 2e[,,)da'da'' 

ik 

where is 

1 / du^ du'^ ^ du' du' ^ 



"ik 



~^ a^i ~^ X! a^i a^k I 



2 V da'' 9a* , 9a* da'' 



The quantity e-jr. is known as the strain tensor. It is fundamental in the theory of elas- 
ticity. In most treatments of elasticity it is assumed that the displacements as well as 
their derivatives are infinitesimal so the last term in Equation ([6]) is dropped. This is an 
approximation that we will not make in this work. 
The quantity 

(9m* du'' ^ du' du' 
= ^-'=+9^+9^ + ?9^9^ 

is the metric for our system and determines the distance between any two points. 
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That this metric is simply the result of a coordinate transformation from the flat space 
metric can be seen by writing the metric in the form^^ 



1 dx^ 


dx' 


dx^ \ 


da^ 


da^ 


da^ 


dx^ 


dx' 


dx^ 


da' 


da' 


da' 


dx^ 


dx' 


dx^ 


\ da^ 


da^ 


da^ I 



( 



1 



where 



1 

Vo ly 

J'^IJ 

du^ 



/ dx^ 


dx^ 


dx^ \ 


da^ 


da' 


da^ 


dx' 


dx' 


dx' 


da^ 


da' 


da^ 


dx^ 


dx^ 


dx^ 


\ da^ 


da' 


da^ / 



da'' da'^' 

and J is the Jacobian of the transformation. Later in section IVI Al we will show that the 
metric for the Fourier modes of our system is not a simple coordinate transformation. 

The inverse matrix ((7*^) = {gik)^^ can be obtained by direct inversion of Equation ([71) 
which would yield components of (?*^ in terms of derivatives of with respect to the internal 
coordinates a*. Later it will be useful however to have the inverse metric stated in terms of 
derivatives of with respect to the external coordinates x*. 

To accomplish this we write the inverse metric as {g"^^) = (J^^)(J^^)"^ where 

/ da^ da^ da^ \ 



dx^ 


dx' 


dx^ 


da' 


da' 


da' 


dx^ 


dx' 


dx^ 


da^ 


da^ 


da^ 



\ dx^ dx' dx^ / 



This yields for the inverse metric 



5. 



du" 



ik 



dx'^ 
5ik + 2eifc 



du'^ _l_ ^ ^""^ ^""^ 



dx 



dx''- dx'^ 



(9) 



and ejfc is the strain tensor in fixed space coordinates. We see that the metric components 
can be written in terms of either sets of coordinates, internal or fixed space. 
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III. EQUATIONS OF MOTION 



In the following we will use the notation 



u 



and therefore the strain tensor is 



1 



We work in the fixed space coordinates and take the strain energy as the lagrangian 
density of our system. This approach leads to the usual equations of equilibrium in elasticity 
theory^ii^. The strain energy is quadratic in the strain tensor e'^'^ and can be written as 



fiuap 



The quantities Cf^i^ap are known as the elastic stiffness constants of the material ^. For an 
isotropic space most of the coefficients are zero and in 3 dimensions, the lagrangian density 
reduces to 



L = (A + 2fi) el^ + + £33 + 2A [eue22 + £11633 + e22e33] + 4/i 

where A and fi are known as Lame constants^. 
The usual Lagrange equations, 



2,2,2 

^12 + ^13 + ^23 



(10) 



E 



d ( dL 



dx'^ \ du 



dL 

duP 



0, 



apply with each component of the displacement vector treated as an independent field vari- 
able. Since our Lagrangian contains no terms in the field u^, Lagrange's equations reduce 
to 

which we shall denote 



Using the above form of the Lagrangian one can write 
dL 



du 



u 



pa j 



pu 



-2A(T - 1[i{upy + Uyf) + A 
+ 2/i 



E + ^^2) ^pv - 2 (en + €22) 



E E '^ISaU/SuSpa — 2 ^ ^auUpa 



U 



pv 



(11) 
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where the divergence of the displacement field is o" = (mh + U22 + mss)- 

The first two terms in Equation f|TT]) are first order in the components u^iy while the last 

terms (in square brackets) are second order and higher in m^j,. 

Let us denote by Ep^ the last two bracketed terms in Equation flTT]) . This allows us to 

write 

V-^-2A^-2.i:^K. + «.)+E^-0. (12) 

The quantity V is a vector and every term in Equation f|T2l) transforms as a vector. As 
such the quantity 

^ dx- 

can always be written as the gradient of a scalar plus the curl of a vector 

^ = Va + V X X- 

This allows us to write 

.^^-2A||-2.(vV.|^).|^.(Vx,r (13) 



or 



where 



= 0, (14) 



(p = [(2A + 4/i)a - a] 

and notation (V x x)'' represents the p component of the vector (V x %)• We see therefore 
that there exists a scalar quantity, in the medium that obeys Laplace's equation. Given that 
a is second order in the strain quantities this implies that reduces to the divergence of the 
displacement field, a, in the infinitesimal strain approximation. 



A. Small Strain Approximation 

In the usual treatment of elasticity theory the strain components and their derivatives 
are taken to be infinitesimal. In this approximation, only the first order terms in any 
equation are kept. We will not make this assumption in this paper. Where necessary we 
will adopt a small strain approximation where the strain components are small but not 
infinitesimal. In using this small strain approximation we will always keep our equations of 
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motion to at least order one order of magnitude higher in the strain components than the 
infinitesimal approximation. 



B. Internal Coordinates 

In sections |V] and IVIt we will need to translate the equations of motion from the fixed 
space coordinates to the internal coordinates. For clarity and to adopt a more consistent 
convention, in the remainder of this text we change notation slightly and write the internal 
coordinates not as a* but as x'* and the fixed space coordinates will continue to be unprimed 
and denoted x*. Now using = — x'* we can write 



^ / dx^ du^ \ d 

2-^ ^rr.i 



9x* c?x* / dx'^ 
d 



dx'^ 



dx'^ J (9x* dx'^ 

Equation (fTSll relates derivatives in the fixed space coordinates x* to derivatives in the 
material coordinates x'*. As mentioned earlier, in the standard treatment of elastic solids 
the displacements as well as their derivatives are assumed to be infinitesimal and so the 
second term in Equation (fT5|) is dropped and there is no distinction made between the x* and 
the x'* coordinates. In this paper we will keep the nonlinear terms in Equation f|T5|) when 
changing coordinates. Hence we will make a distinction between the two sets of coordinates 
and this will be pivotal in the derivations to follow. 



C. Vector Transformation 



The partial derivatives in Equation (fT5|) transform as covectors. Therefore from Equa- 
tion f|T5l) we see that any covector transforms as 

^.=^;-e|^^' (16) 
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when going between the unprimed Euchdean space and the primed coordinate system. In a 
similar manner the components of a vector change as 



dx'' = 




dx^ 



= E + 



3 V 



dx^ 



du' 



) 



dx^ 



or for any vector 



yii ^yi ^yi 




(17) 



where the "up" and "down" notation is used to distinguish components of a vector from 
those of a covector. These transformations will be useful in later sections. 

D. Decomposition of Strain 



The quantity u^^, = l/2{u^y ~Uy^) represents a local rigid body motion of the mediunt^^. 
The quantity e^y = \/2{u^^ + Uy^) represents what is usually termed "pure deformation" 
and for sufficiently small Ufj,u we have e^y ~ e^i, . This decomposition will be useful later in 
identifying the electromagnetic field vector. 

We will now demonstrate a new method for reducing Laplace's Equation ( |T4l) to Dirac's 
equation and compare this method to the traditional Dirac reduction. 



The transformation in Equation (1171) has the form V ^ V — 6V where 




This can be decomposed as^: 




(18) 
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IV. CARTAN'S SPINORS 



The concept of Spinors was introduced by Eli Cartan in 1913^^. In Cartan's original 
formulation spinors were motivated by studying isotropic vectors which are vectors of zero 
length. In three dimensions the equation of an isotropic vector is 

{x^f + {x^f + {x^f = Q (19) 

for complex quantities x*. A closed form solution to this equation is realized as 

x' = -2ioii. x'=^{eo+e,l and x' = ^ - (20) 
where the two quantities are 



/ x^ — ix"^ , ^ / —x'^ — ix"^ 



^0 = ±\/ 7. and ^1 = ±\ 



2 " V 2 

The two component object = {^o,C,i) is a spinor— and any equation of the form f|T9|) has 
a spinor solution. 

In the following we use the notation = d/dx^ and Laplace's equation is written 

:)2 , q2 , q2 



This equation can be viewed as an isotropic vector in the following way. The components of 
the vector are the partial derivative operators d/dx^ acting on the quantity cf). As long as 
the partial derivatives are restricted to acting on the scalar field it has a spinor solution 
given by 





d 


dx^ J 




. d ^ 


l = - 







and 

^1 = ~2 i^'o^ ^ nl2 ) ^ nlT (22) 

where 

D S • 2 1 1 S • 2 

z = X + IX and z = —x + ix 
and the "hat" notation indicates that the quantities ^ are operators. The equations 

so 



dz^ 
and 

~ dz^ 
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are equations of fractional derivatives of order 1/2 denoted .^o = D^Ja^ and = D^Ji^ . Frac- 
tional derivatives have the property thalr^ 

' ' dz 

and various methods exist for writing closed form solutions for these operators^^. The exact 
form for these fractional derivatives however, is not important here. The important thing 
to note is that a solution to Laplace's equation can be written in terms of spinors which are 
fractional derivatives. 

A. Spinor Properties 

If we assume that the fractional derivatives and commute then we also have 

Ml? = ^0^06^1 
- ^ 9 

= -7 (^3 - ^d2) (ds + 182) 
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Using this result combined with Equations fl?I]) and fl22]) we may write for the components 
of our vector 

^ -266 (23) 



dxi 

^ <il + il) (24) 



dx 



2 



Upon complex conjugation we observe that 

{HY = -il and {it)* = -il 

which implies 



^ il-il- (25) 



to = it and ti = iio ■ (26) 
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Under complex conjugation therefore our vector becomes 



\ a I o \ o , I o \ o 

and 



\dx^ ) dx^^ \dx'^ J dx'^ \dx^ J dx^' 

We see therefore that our solution is consistent only for an elastic medium embedded in a 
pseudo-Euclidean space in which the vector d/dx^ is pure imaginary. If we write 

d d 



dx^ ^dt 

then Laplace's equation becomes the wave equation 

We will continue to work with dx rather than dt to avoid carrying the minus sign in the 
computations. 



B. Matrix Form 



It can be readily verified that our spinors satisfy the following equations 

. d 



and in matrix form 



The matrix 



Co 



d 



d 



dx^ dx^ 
I d d 



dx"^ J 
d 

^ dx^ 



dx^ 
d 



— I- 



d \ 



V dx^ ~'~ * dx'^ 



d 



dx^ dx'^ 
d 



dx^ I 

d 



= 

. 9 \ 



dx^ 
d 



dx^ dx"^ 
d d 



dx^ / 



(27) 



V dx^ dx"^ 

is equal to the dot product of the vector = d/dx^ with the pauli spin matrices 



X 



_d 
dx 



d 



d 

dx^ 
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where 



7 




7 



-i 

1 



7 




are the Pauh matrices which satisfy the anticommutation relations 

where / is the identity matrix. 
Equation (!27|) can be written 

E d,Y^ = 0. 



(2^ 



(29) 



where we have used the notation C, = C,4'- This equation has the form of Dirac's equation in 
3 dimensions for a noninteracting massless field ^. 



C. Relation to the Dirac Decomposition 

The fact that Laplace' equation and Dirac's equation are related is not new. However 
the decomposition used here is not the same as that used by Dirac. The usual method of 
connecting the second order Laplace Equation to the first order Dirac equation is to operate 
on Equation fl29l) from the left with J2l=i giving 

= E ^(7V + 7V)5.5.^(x) 

IJ.,U=1 

= {dl + dl + dl)nx) (30) 

where \I' = («i,a;2) is a two component spinor and Equation (!28|) has been used in the last 
step. 

This shows that Dirac's equation does in fact imply Laplace's equation. The important 
thing to note about Equation fl5U]) however, is that the three dimensional Dirac's equation 
implies not one Laplace equation but two in the sense that each component of the spinor 
satisfies this equation. Explicitly stated, Equation (!30|) reads 



/ 



/ 



df + di + di ai 

1^ dl + dl + dl)[a, 
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for the independent scalars ai, a2- 

Conversely, if one starts with Laplace's equation and tries to recover Dirac's equation, it 
is necessary to start with two independent scalars each independently satisfying Laplace's 
equation. In other words, using the usual methods, it is not possible to take a single scalar 
field that satisfies Laplace's equation and recover Dirac's equation for a two component 
spinor. 

What has been demonstrated in the preceding sections is that starting with only one 
scalar quantity satisfying Laplace's equation Dirac's equation for a two component spinor 
may be derived. Furthermore any medium (such as an elastic solid) that has a single scalar 
that satisfies Laplace's equation must have a spinor that satisfies Dirac's equation and such 
a derivation necessitates the use of fractional derivatives. 



V. TRANSFORMATION TO INTERNAL COORDINATES 

In section IVI Al we will take the x'^ coordinate to be periodic and we will derive equations 
for the Fourier components of our fields. Since the elastic solid is assumed to be periodic 
in the internal coordinates we need to translate our equations of motion from fixed space 
coordinates to internal coordinates. Using Equation (fT5|) we can rewrite Equation (129|) . as 

Er(^,-E'£K)(-o (31) 

or 

where d'^ = d/dx'^ and 'j'^ is given by 

7"^ = 7'^-E«M«7"- (32) 

a=l 

By Equation (fTTl) these are the gamma matrices expressed in the primed coordinate system. 
The anticommutator of these matrices is 

{7"^, r} = {7^ - E ^M«7", 7^^ - E «^/57''} 

a /3 

= {7^ Y} - E «^/3{7^ 7^} - E «M-{7", r} + E ^,o.u.p{r, 7^} 

13 Oi al3 
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P a a f3 



2/ I Sf,^ - -u^y + Y 



= 2Ig 

This shows that the gamma matrices have the form of the usual dirac's matrices in a 
curved spaced. To further develop the form of Equation fl3Tl) we have to transform the spinor 
properties of ^. As currently written ^ is a spinor with respect to the Xj coordinates not the 
x[ coordinates. To transform its spinor properties we use a real similarity transformation^ 
and write ^ = where S* is a similarity transformation that takes our spinor in to a 
spinor in x'^. 

We then have 

Equation (13T|) then becomes 

= I'nSd'.^+id'^S)^] 

= s-ws[d'^^ + s-\d'^s)^] 

Using {d'^S~^)S = —S~^{d'^S). This can finally be written 

rK-r,]i = o (33) 

where = {d'^S-^)S and Y = S-^S. 

Equation (l33l) has the form of the Einstein-Dirac equation in 3 dimensions for a free 
particle of zero mass. The quantity d'^ — is the covariant derivative for an object with 
spin in a curved spaced. In order to make this identification, the field must satisfy the 
additional equation^i^ 

where F^g*^ is the usual Christoffel symbol. That this equation holds can be seen by considering 
the equation 

dY _ Q 
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true in the unprimed coordinate system. But since the unprimed coordinate system is 
Euchdean space, the Christoffel symbols are identically zero. This allows us to write 

Since this is a tensor equation true in all frames, in the primed coordinate system we can 
immediately write 

dli>^ + i^Tl = 

Using 7'^ = S'7^5'~\ we have 

d'^isrs-') + {sfs-')rl = 

or 

idis)rs-^ + sidir)s-^ + sriKS'^) + is^^s-^)T'^, = 0. 

Multiplying by S^^ on the left and S on the right yields 

s-\dis)r + {dir) + ndis'')s + = 

Finally, using Ti, = {dlS~^)S and again noting that d'^S'^S = —S'^d^S we have, 

rr, - + {d'.r + 1%'.) = (34) 

We have just demonstrated that in the internal coordinates, the equations of motion 
of an elastic medium have the same form as the free-field Einstein Dirac equation for a 
massless particle in three dimensions. Of course this form is trivial in the sense that it is due 
solely to a coordinate change and can be removed by simply changing back to the unprimed 
coordinates. In the next section we will show that if one of the dimensions in our problem 
is periodic, then the equations of motion for the fourier modes is not trivial. Furthermore 
the introduction of the fourier components will generate extra terms in Equation (I3T]) that 
implies a series of equations relevant for particles with mass coupled to fields that can be 
associated with electromagnetism and classical gravity. 



VI. INTERACTING PARTICLES WITH MASS 



In this section we again consider a three dimensional elastic solid but we take the third 
dimension to be compact with the topology of a circle. All field variables then become 

16 



periodic functions of x'^ and can be Fourier transformed. The act of Fourier transforming 
the equations of motion will effectively reduce the dimensionality of our problem from three 
dimensions to two dimensions This dimensional reduction will in some 

cases create two of the same type of objects (for instance a two dimensional metric and a three 
dimensional metric). We will need to take care to distinguish between the two dimensional 
and three dimensional quantities and use explicit limits in most summations. 

Throughout the derivation we will assume the small strain approximation of section 
(IIII A|) . In doing so, we will keep the equations of motion of the system to one order of 



magnitude higher in the strain components than the infinitesimal results. For the Dirac 
equation the infinitesimal result would yield 'y^d^'^ = which is a linear operator that is 
zero'th order in u^,^, acting on In our derivation extra terms will appear in this equation 
of motion and we will keep the linear operator to first order in Uf^^. 

A. Fourier Transform 

In preparation for Fourier Transforming we isolate the terms involving x'^ and rewrite 
Equation (|3T|) as, 



Similar to Equation fl32l) we define 



7'^ = 7^ - E «M/37'' 

/3=1 



to obtain 



E 7'^ - ns,&^ i - ( E l'u,pus,&\ i + l'U-Y^ ^.3^: - uss&\ i = 

M=l \/3=l J \ u=l J 

In keeping with the small strain approximation the second term in parenthesis may be 
neglected as being second order in the small strain quantities. We are left with 

E 1" - ^3.5^) e + 7' f 5^ - E «^35: - u,,d'\ i = 

Next we transform the spinor 
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Which yields 



E {s [a; - u,,d',] + [{&^s) - u,,{d',s)\) I 



In the absence of deformation there is no distinction between primed and unprimed coor- 
dinates and the similarity transformation therefore reduces to the identity matrix as the 
displacement field Un goes to zero. This imphes that S has the form 

S = I + f{u^) 

where f{u^) is a matrix that is a function of the strain components and / — as the strain 
components go to zero. This gives 



Uapd'^S = Uaj3- 



dx 



a/3 



df duu 
^ du^dx'^ 



E 



which implies that these terms are second order in the strain components and can be ne- 
glected. We are left with 

2 



T.i"s{&,-u,,&, + s-\d'^s))l 

+7''^ [d', E u^^d', - U33% + S-\%S)^ I = 
Now we multiply by —iS~^^^ on the left to obtain 



^3 - E ^^^91 - ussd', + S-\d',S) ]i=Q 



(35) 



where 



i/=i 



r 



This translation to the internal coordinates has introduced several fields into the equations 
of motion. In addition to the spinor field ^ there are the fields u^i, and S~^{d'ij^S). We are 
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now in a position to Fourier transform these field quantities. We first transform tlie fields 
u^y and S^^{d^S) in equation fl35|) to obtain 



u 



and 



s-\d,s) = [s-\d,s)\ 

where u,yf^^k and [5*^^ (9^5*)]^ are the k^^ Fourier modes of the relevant quantities and k = 
27m/ a with a the length of the circle formed by the elastic solid in the Xg direction and n is 
an integer. Equation fl35|) now becomes, 



1^=1 



Next we transform the spinor (remembering that it is periodic in An) 



and the gamma matrices 



r = E7re' 



with q = 2TTj/a, I = 2TTj'/a and j,j' integers. This yields, 



q,k,l 



E 1^ ( - i^u,,^k + s-\d'^s) ^ ] e,/2 







This equation is true for each distinct value of I + k + q/2 = m/2 with m = 2T{n/a and n 
an integer. Writing q = {m — 2k — 21) yields finally. 



E 

Ik 



^~n(^, {m-2k-2l) 

E ^Afi ~ ^ n ^3/.,fc + 



S (d' S) ^(m-2k-2l)/2 



^ ({m-2k-2l) ^ 2 (m-2A;-2/) 

+ dlfi I ^ 4,0 - « Wi/3,fcCi. ^ M33,fc - « 



This is an infinite series of equations describing the dynamics of the fields 
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B. Spectrum of Lowest modes 



In quantum mechanical systems a low energy approximation can often be made in which 
only the lowest Fourier modes of a system are present at low energies (or low temperatures). 
For instance, in condensed matter systems the mean number of phonons of wavevector k 
present in a system at temperature T is 

1 



(37) 



where h is Planck's constant, ks is Boltzmann's constant, T is temperature and Uk for 
acoustic phonons can be approximated by Uk = ck with c a constant. The important aspect 
of Equation (1371) is that as T approaches only the k = mode is occupied. So for spin zero 
fields (such as u^u,o) in the low energy limit in a quantum system only the k = modes are 
present. In this paper we are not considering a quantum mechanical system. Nevertheless 
we shall examine this "low energy" approximation but will not attempt to give a rigorous 
justification, focusing instead on examining the form of the equations of motion that result. 

We therefore consider a theory in which only the lowest k and / modes are present and 
Equation (jSGj) reduces to 

2 / 2 



2 



u=l 

m 



m 



true for each value of m. 

We are now in a position to examine in detail the form of the equations for the spinor 
fields irn/2- 



1. mi mode 

We use the notation niq = 27rg/a corresponding to the equation of motion of the fields 
For clarity we focus on the equation of motion of the q = 1 field which yields 



E 7o" - ^^^3.,0 + [5-^(5;^)] J 6/2 

+ + ^ E «.3,o9: - ^«33,o - ^ [S-'id'.S)]^"^ 6/2 = 



(3^ 
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Let us focus on the first three terms in Equation (l38l) 



(39) 



This portion of the equation of motion is very suggestive when compared with Equation ([T]) 
and suggests we associate with the 2 dimensional gamma matrices and that we associate 
Ms^^o and S~^{d'^S) ^ with the electromagnetic vector potential and the spin connection. In 
order to make these identifications we must show that they are consistent with Maxwell's 
equations, and the auxiliary Equation ([3]). We will now demonstrate that these definitions 
are consistent. 



C. Spin Connection 



We propose to make the following definition for the 2D spin connection 



(40) 



where u;3^ = (1/2)('U3^— -u^a). It must now be shown that this is consistent with Equation (jH]). 
A solution to Equation ([3]) is unique only up to an additive multiple of the unit matrix^^^°. It 
is sufficient therefore to consider only the first term on the right hand side of Equation fHOj) . 
First let us note that the gamma matrices satisfy the anticommutation relations 

Upon Fourier transforming the gamma matrices this becomes 

k,k' k,k' \ (3=1 

and using the ansatz that only the k,k' = modes are present yields 

2 



{70 , 7o } = 2/ I - (u^ufi - u^^fi) + E u^pfiUupfl I . 

13=1 



(41) 



If we insist that the metric of our system is equal to twice the anticommutator of the gamma 
matrices then we are led to define 



21 g^-^{Y^^-] 



(42) 
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which becomes 

21 g"-" - {7o",7o} 

2 

= V ^ i^^^'^,o + '^ufifi) + '^^^^B,oUuf3fi■ (43) 

/3=1 

This is the metric for our two dimensional subspace and it does not have the form of a 
simple coordinate transformation on a flat space metric like that of section [III In other 
words there is no transformation involving the coordinates x'^ and x'^ that will remove the 
Fourier transform in this equation. 

The auxiliary Equation ([3]) can now be shown to hold in a similar manner to the construc- 
tion in three dimensions. First we note that the definition of the two dimensional metric 
Equation (I43p can be obtained from a transformation on 2D Euclidean space 



/3=1 



/ dx^ dx^ \ 

dx'^ dx'^ 

dx^ dx"^ 

^ dx''^ dx''^ ' 



J'^IJ 



I dx^ dx^ \ 

dx'^ dx'"^ 

dx'^ dx"^ 

^ dx'^ dx'"^ ' 



(44) 



(45) 



where 



dx^ 



(46) 



dx"^ ^ dx'' 

In other words the two dimensional metric in Equation (143!) is the Fourier Transform of 
a matrix obtained by a coordinate transformation from two dimensional Euclidean space. 
Given this transformation we note, similar to Equation (fT71) . the components of any vector 
V'^ transform in the same manner as the differential element dx'^ or 

13=1 

We now start with the following tensor equation valid in the Euclidean frame 



dx'' 



/3=1 
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which imphes 



13=1 



or 



/3=1 

vahd in the primed coordinate system. We now use —27^7'^ = S'^^S^^ to obtain, 

-{d',s-^)sr + 7'^(5:^-')5 + i^.i' + E 7^rjt j = 

We now Fourier transform the fields u^i, and {dyS^^)S and use the ansatz that only the 
lowest fourier mode is retained. This gives 



{d,S-^)S ^ \{d,S~^)S 



7o 



(3=1 

This demonstrates that Equation is satisfied. 

We next examine Maxwell's equations and the Einstein Field equations by comparing our 
dimensional reduction to that used in Kaluza Klein theory. 



VII. KALUZA KLEIN THEORY 

In Kaluza Klein Theories one starts with an dimensional system (usually 5 dimen- 
sional) and uses dimensional reduction to reduce the dimensionality of the system to — 1 
dimensions. The starting point is the assumption that the Ricci tensor is identically zero in 

dimensions where the Ricci tensor is defined as 

(jrP f)rp 

p _ ap ap I per J^p J^a pP (A^\ 

~ dxP dx(^ ^ap^Pa K"^' ) 

and are the Christoffel symbols. 

The ansatz is then made that all field variables are independent of the A^th coordinate 
effectively reducing the dimensionality of the problem to A^ — 1 dimensions. The key outcome 
of this dimensional reduction is that the equation Raf3 = in A^ dimensions now becomes 
three sets of equations in A^ — 1 dimensions. Two of these three equations has the form of 
Maxwell's equations and the Einstein Field equations. 
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In the next section we will compare our dimensional reduction with the Kaluza Klein 
methods. In doing so we will continue to use the small strain approximation of Section fllll A|) . 
Our approach will be to keep the field equations to at least second order in the strain 
quantities. From Equation fl47l) this implies that we need to keep the Christoffel symbols to 
second order in Ufj_^. Given the form of the Christoffel symbols 

^ [ dxf^ dx^ dx^ j 
we see that it is sufficient to keep the metric tensor g^p to second order in and the inverse 
metric g""^ to first order in Uap- 



VIII. THE EINSTEIN FIELD AND MAXWELL'S EQUATIONS 

In preparation for a comparison with Kaluza Klein theory, we wish to invert the matrices 
defined in Equation ^ and Equation (11511 keeping only terms that are to second order and 
lower in m^,^. 

In the following we will use the notation that 3 dimensional quantities are denoted with 
a carat (e.g. g^^j) and two dimensional quantities do not have a carat. Our two dimensional 
inverse metric from Equation (j43l) is 



where 



and /i, V take values 1 or 2. 

The inverse of Equation fHHl) is 

9i,v -- 



'l + 2en -2ei2 ^ 



2eio 1 + 2e 



22 / 



1 

D 



( 1 + 2e22 -2ei2 \ 



y -2ei2 1 + 2eii ) 
The determinant of the inverse metric g^^ is 

D = l + 2(en + £22) + 4(eiie22 - e^) 



(48) 
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Using the expansion 



1 - X + + 0{x^) 



we have 



- = 1 - 2(en + €22) + 4(ene22 + 612 + e?i + 4) + 0{€^) 



This gives the following form for the metric tensor correct to second order in u^v 

2 



(49) 



a=l 



In a similar manner the three dimensional metric in Equation ^ can be inverted to yield 



(50) 



with 



a=l 



fi=\ 



and /i, z/ taking on the values 1 to 3. The combination of Equations ([2]), (HUj) and (15UI) 
allows us to write the three dimensional metric in terms of the two dimensional metric 



Qafi + 4630-63/3 — 2eQ,3 + 4 ^vofi-vZ 



9ap 



-2e>3 + 4 ^ e^/3e^3 1 - 2e33 + 4 ^ e 



2 

i/3 



+ 0(m=^ 



and 



9 



a/3 



-Ma3 - M3a + 51 



3 



-M/33 - M3/3 + XI ^yP^y^ 1 ~ 2m33 + X 



1/3 



+ 0(U=') (51) 



we can now compare this dimensional reduction to the usual Kaluza Klein result. The 
Kaluza Klein dimensional reduction from dimensions to — 1 dimensions takes the form 



f 9ap - "^^A^ -<I>M„ \ 



g{kk) 



a/3 



-A' 



\ -A^ (-$2 + AM J / 



(52) 
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where A is the electromagnetic vector potential and $ is a scalar. Comparison of Equa- 
tions (ISTl) and (l52l) suggests that we make the following definitions 



13=1 



$^ = - |l-2e33+4^6j3 

/3=1 



which implies 



^aS^a^ I 

fi=l \ 13=1 J \ a=l / 

2 



(53) 
(54) 



(55) 
(56) 



and 



- $ 



1 + (-2e33 + 4E^=ie 



l + 2e33-4Xe 



2 

/93 



/3=1 



This now allows us to put Equation flSTj) into the same form as the Kaluza Klein decom- 
position, 



9al3 



(57) 



-al3 



+ 0{u'). 



{hi 



\ -A^ (-$-2 + A^A^) j 
The final form is obtained by taking the Fourier transform and keeping only the zero'th term 
of u. 



A" I U3a,0 + Ua3,0 " X "/3a,0M/33,0 

P=l 
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2 

/3=1 

^ 2^33,0 +4E 6^3,0 

with 

^ V /3=1 

The Christoffel symbols and the Ricci scalar can now be written in terms of the two di- 
mensional metric, the electromagnetic vector potential and the scalar field 0. The details 
of this reduction can be found in reference^^ and here we simply quote the results. Using 
Equation ( !57l) . the Christoffel symbols are 



2 

1 



r^3 = ^<^'F^ + 

.-1 



+ A^<l>„) (59) 



and the components of the Ricci tensor are, 

i?44 = ^^'F^/'F^fS + $<l>f. 



where we have used the notation where a comma indicates an ordinary derivative and a 
semicolon indicates covariant differentiation. 

Using the fact that Ra/3 = gives now three sets equations for the scalar (p, the elec- 
tromagnetic field A^ and the 2D Ricci tensor Ra/3- We again quote the results of this 
transformation^^ . 

1 

4 
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These equation can be interpreted as a wave equation for the quantity $, Maxwell's 
equations for the fields and the Einstein field equations, the latter of which may be 
written 



where the quantities T^Jl and T"^ are effective energy momentum tensors given by 



These results indicate that the definition of the electromagnetic vector potential given in 
Equations is consistent. 

IX. DIRAC'S EQUATION 

Given the definitions of the spin connection fl4UI) and the electromagnetic field vector fl53l) 
we now pick up the derivation of Dirac's equation from Equation fl38l) which reads 




(61) 





(62) 



Using Equation (fT8|) the quantity can be written 



and the quantity is 



= -A^-u;3^ + 0{u^) 
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We also write the spinor field as 




Equation fl62|) now becomes 



(63) 




(l-M33,o) + -^^3^<9; 




A. 



Mass Term 



The term ^(1 — Msa^o) can be written in terms of the scalar $^ 



as 



^(l-U,3,0) = ^(l-<l.VO(«') 



This scalar term can be interpreted as a mass term in Equation fl63|) for a particle with mass 
m, where m is 



With the exception of the last three terms, this equation has the same form as Dirac's 
Equation ([T]). This is the main result of this work which shows that a Dirac-like equation is 
present in this system. 

B. Additional Interaction Terms 

The last three terms in Equation fl64p can be interpreted as interaction terms between 
the spinor field \1/ and the vector fields A", and u^u, and the matrix field [S'~^(93S')]q. These 
terms do not normally appear in Dirac's equation but the existence of these types of terms 
is not forbidden in modern quantum field theories^. 




Equation (jHS]) now becomes 




(64) 
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In addition to the equation given in (16^ . each of these fields satisfies additional constraint 
equations that determines its dynamics. The electromagnetic field satisfies Maxwell's Equa- 
tion and has been discussed in section IVllII The interaction with the matrix [S'^^(93S')]q 
results from the spin connection in 3 dimensions and its solution is determined from Equa- 
tion (1341) ■ An equation of motion for the field uo^^^ can be derived from Equation (IT^ . If we 
take the divergence of the vector from this equation we have 



which implies 



VV=-(2A + 2,)|^ + |i + (Vx,r 



We can now write the equation of motion of the field ujpp as 

- (15 - 3 = ^■>' - ^^"^ 

We see that the scalar (and hence the spinor field ^) is completely determined by the 
scalars a and a in Equation (fT51) and the field Upi, is completely determined by the vector 
field X- 

We will not attempt to identify these three additional terms with known interactions, 
simply noting again that interactions of this type are not forbidden to exist and that the 
central result of this work is Equation flMl) which shows that there is a Dirac-like equation 
that exists in a classical elastic solid. 



X. HIGHER DIMENSIONS 



Throughout our derivation we assumed that we were working in three dimensional space. 
This was convenient since explicit solutions to Laplace's equation could be obtained in terms 
of fractional derivatives. The formalism however extends to any number of dimensions^^ii^. 
The main difference is that there is no explicit solution of the spinors in terms of fractional 
derivatives. Nevertheless it should be noted that the basic equations hold in higher dimen- 
sions and therefore the usual four dimensional Dirac equations could in principle be derived 
by starting with an elastic solid in 5 dimensions. 
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XI. CONCLUSIONS 



We have taken a model of an elastic medium and derived an equation of motion that 
has the same form as Dirac's equation in the presence of electromagnetism and gravity. 
We derived our equation by using the formalism of Cartan to reduce the quadratic form 
of Laplace's equation to the linear form of Dirac's equation. We further assumed that 
one coordinate was compact and upon Fourier transforming this coordinate we obtained 
a mass term and an electromagnetic interaction term in the equations of motion. The 
formalism demonstrates that the classical Einstein-Dirac-Maxwell equations can be derived 
as the Equations of motion of the Fourier modes of an elastic sohd in the small strain 
approximation. 
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